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Quantitative features of multifractal subtleties in time series 
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PACS 05 . 45 . Tp - Time series analysis 

PACS 05 . 40 . Fb - Random walks and Levy flights 

PACS 05 . 45 . Df - Fractals 
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systems 

Abstract. - Based on the Multifractal Detrended Fluctuation Analysis (MFDFA) and on the 
Wavelet Transform Modulus Maxima (WTMM) methods we investigate the origin of multifractal- 
ity in the time series. Series fluctuating according to a (/Gaussian distribution, both uncorrelated 
and correlated in time, are used. For the uncorrelated series at the border (5 = 5/3) between the 
Gaussian and the Levy basins of attraction asymptotically we find a phase-like transition between 
monofractal and bifractal characteristics. This indicates that these may solely be the specific 
nonlinear temporal correlations that organize the series into a genuine multifractal hierarchy. For 
analyzing various features of multifractality due to such correlations, we use the model series gen- 
erated from the binomial cascade as well as empirical series. Then, within the temporal ranges 
of well developed power-law correlations we find a fast convergence in all multifractal measures. 
Besides of its practical significance this fact may reflect another manifestation of a conjectured 
g-generalized Central Limit Theorem. 
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Multiscaling represents a concept that indicates an ex- 
tremely promising direction towards grasping the essential 
characteristics of complexity. Indeed, the related measure 
in terms of multifractal spectra offers an attractively com- 
pact frame to quantify the hierarchy of scales and speci- 
ficity of their interwoven organization. This in particular 
applies to the dynamical aspects of complexity. The em- 
pirics of such aspects is typically accessible through the 
time series containing information about the consecutive 
, states of a complex system. A plenty of evidence has been 
reported that empirical data coming from such diverse 
fields as physics of turbulence [I], geophysics [2], astro- 
physics [3] , physics of plasma |4] , physiology [5] , complex 
networks research [B] and econophysics [7] present signa- 
tures of multiscaling. As far as the time series are con- 
cerned an issue that however still remains open in the 
literature is which identifiable distinct attribute of com- 
plexity is in this context relevant. It seems commonly 
believed that either of the two, the nonlinear temporal 
correlations as well as abundantly accompanying them the 
non-Gaussian heavy tails of fluctuations, or both equally, 
may account for the multiscaling effects. Bearing in mind 
the very definition of the multifractal spectrum [9] through 



its relation to the strength and organization of singulari- 
ties in a series and a more recent work [TU] devoted to the 
construction and analysis of continuous multifractal cas- 
cades it seems however more natural to expect that these 
are the temporal correlations that are most relevant and 
even their scale free character demanded. 

In order to quantitatively illuminate on the related is- 
sues we impose a selective discrimination by elaborating 
first of all the series with no time correlations at all, but 
instead with the systematically varying degree of heavi- 
ness of tails in their fluctuations. As an economic and 
physically motivated choice we use series sampled from 
the ^Gaussian distribution 
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whose origin can be routed to the nonextensive generaliza- 
tion of entropy [TTlfTS] . The special q = 1 case corresponds 
to the standard Gaussian distribution. 

Further reason for considering this family of distribu- 
tions is that they often well model the empirical data 
with heavy-tailed pdf's [T21Q3]. Moreover, according to 
the generalized Central Limit Theorem (CLT), the uncor- 
related ^Gaussian signals for q > 1 are unstable under 
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convolution, and, depending on q, they are either in the 
standard Gaussian (1 < q < 5/3) or in the Levy stable 
distributions (5/3 < q < 3) basin of attraction with re- 
lation a L — (3 — q)/(q — 1) [II]- From this perspective 
an anticipated asymptotic for the uncorrelated series with 
(/Gaussian pdf is a monofractal (Gaussian attractor) or a 
multifractal (Levy attractor) [13] which for the finite se- 
ries assumes a bifractal form as for the truncated Levy 
processes [IB] . 

At present the Multifractal Detrended Fluctuation 
Analysis (MFDFA) [17] and the Wavelet Transform Modu- 
lus Maxima (WTMM) [18] are the two most efficient prac- 
tical methods to quantify multifractality. Both have their 
advantages, but the comparative evaluation of their over- 
all performance [19] indicates that in most cases the more 
reliable choice is MFDFA. Therefore we use this method 
as our basic tool and occasionally WTMM (according to 
the procedure described for instance in ref. |19| ) as an 
auxiliary independent consistency test. 

The main points of MFDFA go as follows. Let 
jv be a discrete signal. In MFDFA, one starts 

with the signal profile Y(j) = ( x (^)~ < x >)> 3 = 

1, N, where < ... > denotes averaging over all i's. Then 
one divides the Y(J) into M n non-overlapping segments of 
length n (n < N) starting from both the beginning and 
the end of the signal (2M n segments total). For each seg- 
ment a local trend is estimated by fitting an Ith order 
polynomial Pj , which is then subtracted from the sig- 
nal profile. For the so-detrended signal a local variance 
F 2 (v,n) in each segment v is calculated for each, from 
n m in to n max , scale variable n. Here, if not stated oth- 
erwise, we take n m i n = 40 and n max — N/5. Finally, by 
averaging F 2 (v : n) over all segments v one calculates the 
rth order fluctuation function: 

f 1 2A/„ s l/ r 

k V— 1 J 

where r E R. Scaling behavior of the fluctuation func- 
tion F r (n) ~ n h ^ is an indication that the analyzed 
signal has fractal structure. If h(r) = const, the sig- 
nal is monofractal; it is multifractal otherwise. In both 
cases the singularity spectrum f(a) [9] can be calculated: 
f(a) — r[a — h(r)] + 1, where a = h(r) + rh'(r) is the 
singularity strength. 

In the WTMM one calculates the wavelet transform 
T^(n, s) of a time series x(i), where ip is the wavelet 
shifted by n, and s is scale. For each scale one finds 
the local maxima of and then evaluates the parti- 
tion function: Z(r, s) = YlieL(s) \T-iii(ni(s), s)\ r , where 
L(s) denotes the set of all maxima for the scale s 
and ni(s) stands for the position of a particular max- 
imum. To preserve monotonicity of Z(r, s) for differ- 
ent values of s, an additional criterion must be applied: 
Z(r,s) = J2ieL(,s)( sn Ps'< s \ T -<p( n i(s'),s')\) r . For fractal 
signals Z(r,s) ~ s T ( r \ The scaling exponent r(r) is di- 



rectly related with a and f{a) by the Legendre transform: 
a = r'(r) and f(a) — ra — r(r). 

Typically, the f{a) spectrum of a multifractal set resem- 
bles an inverted parabola which for monofractal signals re- 
duces to a single point. In practice, however, one usually 
obtains a very narrow parabola even for the monofractal 
signals. The degree of multifractality can be evaluated by 
measuring width |17) 

Aa = a max - a min (3) 

between the calculated extremal values of a for a given 
range r min < r < r max . 

In order to assess the impact of a finite signal length on 
/(a) for the gGaussian uncorrelated series, we generated 5 
independent realizations of such series of different lengths: 
N = 10 4 ; 10 5 ; 10 6 , for a sequence of q values from the range 
1 < q < 2 with a step Aq = 0.05. This interval thus covers 
also the crossover regime near q = 5/3 between the Gaus- 
sian and the Levy attractor. The second- order polynomial 
as the detrending filter in MFDFA and the third deriva- 
tive of Gaussian as the mother wavelet in WTMM are 
used. As we consider signals with heavy tails, the range 
of the index r cannot be too broad; we choose r 6 [—4, 4]. 
This relatively narrow range of r values allowed, automat- 
ically prevents a potential distortion of the results by the 
so-called " freezing" phenomenon which manifests itself in 
the linearization of h(r) at large r as observed for a log- 
normal cascade 21,22 . This linearization may thus lead 
to an artificial narrowing of the singularity spectrum. In 
this connection it is also worth to notice - as a side remark 
- that as ref. [19 indicates the MFDFA procedure is less 
susceptible to such effects than WTMM. 

F r (n) (MFDFA) calculated for a few exemplary choices 
of q between 1.5 to 1.8 are shown in Figure 1 for N = 10 6 . 
For q = 1.50 (top left) there is a good power-law depen- 
dence of F r (n) for the whole range of scales and a small 
variability of the slope coefficients indicating a monofrac- 
tal character of the analyzed data. On the opposite edge 
(q > 1.70, bottom row) one also observes a good scal- 
ing for all the scales and all the values of r, but here 
the slope coefficients are r-dependent. Between these two 
cases (1.55 < q < 1.65) there are signals with no uniform 
scaling of F r (n) . Instead, we can distinguish two regions 
of different character: the small- scales region with the 
r-sensitive slope coefficients and the large-scales region 
without significant dependence on r. When going from 
q = 1.55 to q = 1.70, i.e. crossing the border between the 
Gaussian and the Levy regimes at q — 5/3, the multiscal- 
ing region gradually extends towards the larger scales in 
the expense of the monofractal region. For shorter time 
series the crossover is more blurred. For the shortest con- 
sidered signals it is at all difficult to distinguish the two 
regimes. One principal reason for these effects is that the 
scale parameter n in Eq. ([2J needs to be sufficiently large 
so that through the requirements imposed on F 2 (v,n) by 
the CLT the fluctuation function sufficiently accurately 
probes the attractor. 
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Fig. 1: Fluctuation function F r {n) (MFDFA) calculated for the 
signals of length N = 10 6 sampled from different qGaussians. 
In each panel F r (n) for different moments indexed by r £ 
[—4, 4] is displayed. 

Figure 2 shows the width Aa(q) calculated with 
MFDFA (main panel) and with WTMM (inset). If two 
scaling regimes are observed in the F r (n) plots, the corre- 
sponding widths are calculated separately for each regime. 
Analytically calculated f{a) spectra for the finite signals 
in the Levy basin of attraction are bifractal and consist of 
only two points located at (0,0) and (l/aL,l) [1 6 . Nu- 
merically, such bifractal spectra cannot ideally be resolved 
and are represented by the shapes as in Figure 3. Thus 
the width in the main panel of Figure 2 is calculated from 
the continuous spectra for all the q- values and it is due to 
such effects that it assumes values slightly larger than the 
theoretical value Aa (Lovy) = l/a L 1!) for q > 1.8. If we 
approach the crossover at q = 5/3, the deflection of Aa 
from A«' Lcvy ' becomes more significant. In the Gaussian 
basin of attraction all the considered signals doubtlessly 
appear monofractal only for q < 1.3, where the CLT con- 
vergence is fast. However, the signals with N > 10 5 are 
long enough for the CLT to exert its influence even for q 
up to 1.5. 

For 1.5 < q < 5/3 the gGaussians are significantly 
heavy-tailed, which makes their convergence to the Gaus- 
sian attractor very slow. This situation corresponds to 
the coexistence of two scaling regimes in Figure 1. The 
closer to the crossover point we are, the larger scale n is 
needed to obtain the monofractal scaling in F r (n). In this 
region there exists no a priori criterion to choose one scal- 
ing regime and neglect the other, hence for such signals 
there is no unique value of Aa. This problem may be re- 
solved only by considering signals of much larger length 
(e.g., N = 10 9 ), which can sometimes be generated for ar- 
tificial signals, but are typically beyond reach for empirical 
data. The case of our shortest signals (N = 10 4 ) illustrates 
how data with no actual scaling can erronously be consid- 
ered suitable for calculating the scaling exponents. Inset 
of Figure 2 shows the same quantity but obtained from 



Fig. 2: The spectrum width Aa as a function of the gGaussian 
parameter q. Dashed lines for q < 5/3 denote results obtained 
in the two competitive scaling regions of F r (n). The vertical 
dashed lines separate the Gaussian and the Levy basin of at- 
traction. (Main panel) Results obtained with MFDFA. (Inset) 
Results obtained with WTMM. Error bars denote standard de- 
viation of 5 independent realizations of the process. 



the wavelet-based analysis. Due to a restricted applica- 
bility of the WTMM method for shorter signals, only the 
N = 10 5 and N = 10 6 cases are displayed. These results 
are consistent with the ones by MFDFA. 

In Figure 3 we show the /(a) spectra for signals from 
the Levy basin of attraction (1.7 < q < 2.0). They look 
multifractal but the high concentration of data points in 
two small regions near the points (0,0) and (l/ai,,0) in- 
dicates the bifractal structure. Inset in Figure 3 displays 
a max (g). Outcomes of the numerical simulations for differ- 
ent N almost follow the analytical prediction for q < 1.5 
(where a max ~ = 1/2) and for q > 1.8. In the vicin- 

ity of the threshold (1.5 < q < 1.8), the numerical value 
of a max is larger than analytically expected. 

Of an even greater interest are the series that involve 
correlations as well. As one example of such series, whose 
correlations are controllable and the asymptotic true re- 
sult [20] is in addition known explicitly, we elect a signal 
generated by the log-normal cascade at the k — 20 stage 
of cascading. This true result is typically well approxi- 
mated already for k > 12. In the present case the cas- 
cade parameter A = 1.1 (the same notation as in ref. [19] 
is used). The temporal correlations in such a series de- 
cay uniformly on all the scales according to the power-law 
and, as using the cumulative ^Gaussian distribution [T3] 
is shown in Figure 4(c), the distribution of fluctuations at 
this stage (k = 20) of cascading fits well to the gGaussian 
with q = 1.6. As before, the TV = 10 4 ;10 5 ;- 10 6 (pre- 
cisely 2 20 in this last case) long series are extracted from 
the same k = 20 cascade (2 fc data points), they thus obey 
the same distribution of fluctuations, r £ [—4, 4] and the 
scale variable n £ [40, N/5]. For N = 10 4 and N = 10 5 
the result of calculations is averaged over the all nonover- 
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Fig. 3: (Main panel) MFDFA calculated singularity spectra 
f(a) for signals with gGaussian pdf's for a few different choices 
of q from the Levy basin of attraction. Analytical spectra, not 
shown here, are bifractal. (Inset) Inverse of the numerically 
calculated maximum a max , which is directly related to the Levy 
parameter ai. 



lapping segments. The final result is shown in Figure 4(a). 
This result turns out essentially identical for different re- 
alisations of the above specified cascade. The convergence 
with N of /(a) to the correct known result can be seen 
to be encouragingly fast. In the next step, destroying cor- 
relations in these series, by randomly shuffling their data 
points, results in fluctuation functions as in Figure 1 for 
q = 1.6 with their characteristically increasing split in r 
for the decreasing scale variables n. Consequently, extract- 
ing here the scaling exponents from F r (n) for n varying 
within 40 up to say 2 x 10 3 interval - as by necessity would 
be done for a few thousand data points long series - re- 
sults in f(a) shown in Figure 4(b). This outcome, largely 
even stable with increasing N, makes impression of reveal- 
ing multifractality. By running the same procedure for n 
from 2 x 10 3 up to the usual N/5 and increasing N results 
in /(a) that systematically, though slowly, approaches an 
anticipated monofractal character. 

As a final example we elaborate the real empirical data. 
This is the time series representing the high frequency, 
15-second logarithmic returns of the German stock mar- 
ket index (DAX) from the period November 28th, 1997 - 
December 31th, 1999 with the trading hours 8:30 - 17:00 
from the beginning which on September 20th 1999 until 
the end of this period was changed to 9:00 - 17:30. The 
first nonzero return for each day is removed from the anal- 
ysis which results in N = 1, 064, 960 data points long time 
series. The corresponding distribution of fluctuations fits 
also well by the qGaussian with g rj 1.6 which is shown 
in Figure 5(c). Exploring and quantifying the multiscal- 
ing characteristics in the financial time series is one of the 
mainstream issues in econophysics. Indeed, for our DAX 
series we obtain a satisfactory scaling of F r (n) which al- 
lows to evaluate h(r). In order to permit making parallels 
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Fig. 4: (a) Singularity spectra f(a) for log-normal multi- 
plicative cascades of different lengths, (b) f(a) for randomly 
shuffled cascade data, calculated for two different ranges of 
scaling variable: 40 < n < 2,000 ("non-CLT regime") and 
20,000 < n < N/5 ( "CLT regime"), (c) Cumulative distribu- 
tion function for the log- normal cascade (circles) with best- 
fitted gGaussian distribution (dashed line). 



with the cases considered above, we also divide the whole 
series into shorter pieces of length N — 10 5 or N = 10 4 and 
average the results over all the pieces of the same length. 
The scaling variable n € [40, N/5] and the result is shown 
in Figure 5(a). At first glance somewhat unexpectedly (in 
view of the results for the log-normal cascade) one sees 
here no stable convergent result for f(a) in function o( N. 
The span of the so-calculated f(a) systematically shrinks. 
Only a deeper inspection of the functional character of 
nonlinear correlations, as expressed for instance by the 
volatility (return modulus) autocorrelation, shown in Fig- 
ure 6(a), opens perspective to understand this result. As 
expected, such an autocorrelation slowly decays accord- 
ing to the power-law. What however is not so expected 
in view of the present understanding of such effects thus 
far documented in the literature is that autocorrelation 
between events that are separated by more than about 
2 x 10 4 basic units (here 15 seconds) suddenly drops down 
and starts oscillating between the positive and negative 
values with a decreasing amplitude (Figure 6 (b)). Such 
distant events are already not so consistently organized 
into a temporally power-law organized hierarchy relative 
to each other. Taking the scaling variable n larger than 
the power-law correlation range thus introduces admix- 
ture of such a noise- like contribution which explains both 
the shrinkage of /(a) and a shift of its maximum towards 
a = 0.5. That this is a relevant effect can be seen from 
Figure 5(b) where f(a) is calculated for n up to 2 x 10 4 
independently on N. Now both N = 10 5 and N = 10 6 
generate essentially the same result. We also verified that 
even going with n up to only 2 x 10 3 already gives very 
similar result as the one that above is found as stable. 
The above results show that multifractality is a valu- 
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Fig. 5: Singularity spectra /(a) for empirical time series (Ger- 
man DAX 1 min returns) of different lengths, calculated in two 
different ranges of scale parameter: (a) 40 < n < N/5 and (b) 
40 < n < 20, 000. (c) Cumulative distribution function for 
DAX returns (circles) with gGaussian distributions (dashed) 
best fitted separately for negative and positive returns. 



Fig. 6: Autocorrelation function for absolute DAX returns 
shown in log- log scale (a) and in the linear scale (b); note the 
region of negative autocorrelation values. 



cation of several by now available [25] statistical methods 
testing multifractality is recommended. 



able but at the same time subtle concept. In the time se- 
ries the genuine multifractality originates from temporal 
correlations. Encouragingly, within the range of well de- 
veloped power-law temporal correlations, even in the most 
"critical" region of fluctuations just at the border between 
the Gaussian and the Levy basins of attraction, already a 
few thousand points long series suffices to properly evalu- 
ate f(a) using MFDFA. As an example with the financial 
data shows such an analysis needs however necessarily to 
be accompanied with evaluation of the range of those cor- 
relations and should be restricted to the intervals they 
consistently behave according to the power-law. In the re- 
alistic cases such correlations are likely to extend over the 
finite time-lag intervals only. Equivalently, a test of sta- 
bility of the spectrum with the varying length of the series 
may serve the same purpose. Interestingly, much longer 
series are needed to reach convergence for the gGaussian 
timc-uncorrclated scries, especially in the region charac- 
terized by q close to 5/3. In this latter case two results 
are asymptotically possible: either a monofractal singular- 
ity spectrum for q < 5/3 or a bifractal for q > 5/3 which 
means that no genuine multifractality is due to fluctua- 
tions. The observed faster convergence in the presence 
of power-law correlations may indicate another manifesta- 
tion of the conjectured ^-generalized CLT ,23,24 . Finally, 
the effects here elaborated and systematized, when taken 
into account, are expected to be crucially helpful in iden- 
tifying the real multiscaling phenomena, their range and 
significance, and to reliably discern them from what often 
may turn out "apparent" . Once a real multiscaling in this 
way is identified, in order to reach a higher level of con- 
fidence on the related quantitative characteristics, appli- 
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